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We use Monte Carlo simulations to examine the simplest model of an ionic liquid, called the 
restricted primitive model, at a metal surface. We find that at moderately low temperatures the 
capacitance of the metal/ionic liquid interface is so large that the effective thickness of the electro- 
static double-layer is up to 3 times smaller than the ion radius. To interpret these results we suggest 
an approach which is based on the interaction between discrete ions and their image charges in the 
metal surface and which therefore goes beyond the mean-field approximation. When a voltage is 
applied across the interface, the strong image attraction causes counterions to condense onto the 
metal surface to form compact ion-image dipoles. These dipoles repel each other to form a corre- 
lated liquid. When the surface density of these dipoles is low, the insertion of an additional dipole 
does not require much energy. This leads to a large capacitance C that decreases monotonically 
with voltage V , producing a "bell-shaped" curve C{V). We also consider what happens when the 
electrode is made from a semi-metal rather than a perfect metal. In this case, the finite screening 
radius of the electrode shifts the reflection plane for image charges to the interior of the electrode 
and we arrive at a "camel-shaped" C{V). These predictions seem to be in qualitative agreement 
with experiment. 
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I. INTRODUCTION 



Ionic liquids are molten salts made from ions which 
are large enough that their Coulomb interaction is rela- 
tively small, so that they remain in a fluid state at room 
temperature. EssentiaUy, an ionic liquid is a solvent-free 
electrolyte, which means that ionic liquids can be ideally 
suited for applications which require a thin or intensely 
concentrated layer of ionic charge. Ionic liquids are al- 
ready being used for batteries and "supercapacitors" [l| , 
as well as for gating of new electronic materials. It has 
therefore become a subject of great interest to under- 
stand the nature of the interface between an ionic liquid 
and a metallic electrode. 

In its simplest form, the question of how an ionic liquid 
behaves in the vicinity of a charged metal surface seems 
remarkably straightforward. While real-life experiments 
probing the structure of the ionic double-layer can be 
marked by a number of complications [3], the essential 
description is encapsulated in a very simple model: an 
infinite, planar, metallic electrode is placed in contact 
with a semi-infinite volume with uniform dielectric con- 
stant e that contains a total concentration N of mobile 
positively- and negatively- charged hard spheres, each 
with the same diameter a and the same absolute value 
of charge e. Such a model of the ionic liquid is called 
the "restricted primitive model" (RPM). If a voltage V 
is applied between the electrode and the bulk of the ionic 
Uquid, how large is the charge density a of the metal sur- 
face? In other words, what is the capacitance per unit 
area C{V) = da/dV of the interface? 

The answer to this question is well-known in the limit 
of low ion density, large temperature, and low applied 
voltage. In this case the ionic double-layer is well- 
described as a diffuse screening layer with a characteristic 



size equal to the Debye-Hiickel (DH) screening radius 
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Here, fc^T is the thermal energy [Eq. ([l}, and the remain- 
der of this article, uses Gaussian units]. The diffuse layer 
of counter-charge effectively comprises the second half of 
a parallel-plate capacitor of thickness r^H, so that the 
capacitance per unit area is equal to Cdh = s/^T^rDH- 
This result for capacitance is valid as long as the ion 
density is low enough that Na^ <^ 1, the temperature is 
high enough that T ^Tq = e^ /kBea, and the voltage is 
small enough that eV <^ ksT. Under this fairly extreme 
set of assumptions, the ionic double-layer is essentially a 
small perturbation of the bulk density N , so that the ion 
density and electric potential can be described using the 
linearized Poisson-Boltzmann equation. 

More generally, one can characterize the capacitance 
by the effective thickness of the double-layer d* = e/AirC. 
In the DH limit, d* = tdh- In realistic situations, how- 
ever, the characteristic temperature Tq is very large (for 
a = 1 nm and e = 3, Tq ~ 5600 K) and the DH approxi- 
mation fails. One can think that as a result ions become 
more strongly bound to the charged electrode and the 
size of the double-layer shrinks, so that d* decreases and 
the capacitance grows. One may ask, then, how thin 
the double-layer can be, or in other words, how large the 
capacitance can be. The apparent answer to this ques- 
tion goes back to Hclmholtz [3| , who imagined that in an 
extreme case a neutralizing layer of ions could collapse 
completely onto the electrode surface, thereby forming 
the second half of a plane capacitor at a distance equal 
to the ion radius a/2. The result is a double-layer of size 
d* = a/2 and a capacitance per unit area equal to the 
"Helmholtz capacitance" 
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In classical mean- field theories of the electrostatic double- 
layer j^-Qi ^■iid in the recent influential theory of the 
metal/ionic liquid capacitance which accounts for the 
excluded volume among ions JT"!, Ch plays the role of 
a maximum possible capacitance per unit area. Monte 
Carlo simulations |8l-[l0l| seem to confirm this statement. 
However, these and the majority of other simulations 
make the simplification of replacing the metal electrode 
by a uniformly-charged, insulating plane. We argue be- 
low that in this way the essential physics of image charges 
in the metal surface is lost (see also Ref. [ll| , where some 
such simulations are critically analyzed). 

It is the purpose of this article to demonstrate that 
capacitance C > Ch is possible, or in other words, that 
the effective thickness of the double-layer can be smaller 
than the ion radius. Our previous work 112] has demon- 
strated that capacitance C > Ch can occur for highly 
asymmetric ionic liquids (where the cation has a much 
smaller radius than the anion, or vice- versa). Here we 
show that even in the RPM, where cations and anions 
have equal diameter, capacitance significantly larger than 
the Helmholtz value is possible. As we describe below, 
the metallic nature of the electrode — specifically, the 
ability of ions to form image charges in the metal surface 
— plays a key role in the development of large capaci- 
tance. We present the results of Monte Carlo (MC) simu- 
lations of the restricted primitive model of an ionic liquid 
at various temperatures and densities, and we suggest a 
basic theoretical explanation of these results based on the 
weak repulsion between dipoles composed of bound ions 
and their images in the metal surface. 

The remainder of this paper is organized as follows. In 
section II we present our MC results for C{T) at small 
voltage and give them a qualitative explanation. In sec- 
tion HI we explain our MC procedure. Section IV is 
devoted to our analytic theory, including both tempera- 
ture and voltage dependence of the capacitance C{V, T). 
Section V considers the role of the electrode material on 
C{V,T). Section VI discusses analogies between non- 
trivial capacitance phenomena at the metal/ionic liquid 
interface and those in semiconductor devices such as sil- 
icon MOSFETs and gated GaAlAs structures. We con- 
clude in section VII with a summary of our main results. 



II. TEMPERATURE DEPENDENCE OF THE 

CAPACITANCE OF THE METAL/IONIC LIQUID 

INTERFACE 

Fig. [1] shows the zero voltage capacitance C{Q)/Ch, 
as measured by our MC simulations, as a function of 
reduced temperature T* = T/Tq for three different di- 
mensionless ion densities Na^ . The points correspond to 
results from the MC simulation, and solid lines are a fit 
to the form C /Ch = A-{T*)~^'^ , where ^ is a numerical 
constant. The motivation for this (T*)"^/'^ dependence is 
explained in section IV. For all three values of the density 
that we examined, the capacitance at low temperature is 



significantly higher than the Helmholtz value. For prac- 
tical applications, Na^ k, 0.5, and for e = 3 and a = 1 
nm room temperature corresponds to T* « 0.06. Using 
these parameters, the maximum value of capacitance in 
Fig. [U is 3Ci/ « 27/iF/cm2. 
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FIG. 1: The capacitance of the metal/ionic liquid double-layer 
at zero voltage as a function of the dimensionless temperature 
T* , plotted for three values of the dimensionless ion density 
No?' . Symbols represent results from the MC simulation and 
solid lines show a best fit to the form C/Ch = A ■ (T*)"^/^ 
for each density. Error bars are smaller than the symbol size. 

These results should be contrasted with previous sim- 
ulation studies [8;, 9;] of the capacitance of the RPM ionic 
liquid, in which the metallic electrode was replaced by a 
charged insulator with uniform charge density a. These 
studies report a capacitance C(0) that grows with de- 
creasing T* before reaching a peak at T* = T* and then 
coUapsing rapidly at T* < T*. For No? = 0.08 and 0.64, 
T* K, 0.17 and 0.28, respectively. 

The collapse of the capacitance at low temperatures 
T* < T* was interpreted by the authors of Refs. [1, Q as 
the result of strong binding of positive and negative ions 
to form neutral pairs. Such binding leads to an extreme 
sparsity of free charges in the ionic liquid, so that their 
total concentration N j <C N . Substituting N f for iV into 
Eq. ([T]) at r* <C 1, we arrive at a large screening radius 
tbh and therefore much smaller capacitance C(0). These 
arguments are generic and convincing. Why, then, does 
the capacitance in Fig.[T]continue to grow with decreasing 
temperature? 

Here we present a qualitative answer to this question. 
We begin by observing that when the electrode is metal- 
lic, the energy of an ion binding to its image charge, 
— e^/2ea, is exactly half the energy of a bound ion pair, 
—e^/ea. This fact implies that if an ion pair is sepa- 
rated in the bulk and then both ions are brought to the 
metal surface there is no net change in electrostatic en- 
ergy. Thus, even in the absence of applied voltage there 
are plenty of free charges at the metallic surface. This 
allows the double-layer to be very thin and leads to the 
large capacitance shown in Fig. [TJ 

We note that the critical role of image charges for the 



structure of the double-layer has in fact been recognized 
by previous authors |13| . who performed similar simula- 
tions which account explicitly for the electronic polariza- 
tion of the electrode. However, Ref. [1J| explored only 
very low ion density Na'^ — 0.01, where C(0) < Ch- An- 
other paper [IJl studied the capacitance of an ionic liquid 
between two identical metal plates and obtained large ca- 
pacitance C(0) ~ 2Ch, but this study used a much more 
complicated model of the ionic liquid. 

The following section explains our Monte Carlo proce- 
dure. 



III. MONTE CARLO SIMULATION 

In our MC simulations, a canonical ensemble of Ma 
anions and Mc cations is placed in a square prism cell 
of dimensions L x L x L/2 and corresponding volume 
fi = L^ /2. The metallic electrode coincides with one of 
the cell's square faces. Every charge within the cell forms 
an electrostatic image in the electrode surface (z = 0), 
i.e. a charge q = ±e at position {x,y,z) has an image 
charge —q located at {x,y,—z). The total electrostatic 
energy £ of the cell is calculated as 1/2 times the energy 
of a system twice as large composed of the real charges 
and their images, so that 



-, Mt 



(3) 



where Mt = 2{Ma + Mc) is the total number of charges in 
the system (ions plus images), dij is the distance between 
particles i and j, and the two-particle interaction energy 
u{di,j) is 



u{dij) 



< a 



q^qj/edij, dij > a 



(4) 



Here, qi = ±e is the charge of ion i. 

The charge of the electrode is varied by changing the 
number of anions Ma and cations M^ in the system by 
equal and opposite amounts, so that the total number of 
ions Ma + Mc = Nfl remains fixed for a given overall 
density N. The corresponding electronic charge (in the 
form of image charges) in the electrode is Q = e{Ma — 
Mc) and the capacitance dQ/dV can be determined from 
the resulting voltage. We use the system size L = 20a 
everywhere. 

At the beginning of each simulation, positive and neg- 
ative ions are placed within the simulation cell in such 
a way that they do not overlap with each other or with 
the metal surface. The MC program then selects an ion 
at random and attempts to reposition it to a random 
position within a cubic volume of (20)^^ centered on the 
ion's current position. The change in the energy £ asso- 
ciated with this move is then calculated, and the move 
is accepted or rejected based on the standard Metropolis 



algorithm. For one in every 100 attempted moves, the 
MC program chooses the random position from within a 
larger volume (lOa)'^ as a means of overcoming the effects 
of any large, local energy barriers. The simulation cell is 
given periodic boundaries, so that an ion exiting one face 
of the cell re-enters at the opposite face. To ensure ther- 
malization, 2,500 moves per ion are attempted before any 
simulation data is collected. After thermalization, simu- 
lations attempt 2 x 10'' moves per ion, of which 15% - 
50% are accepted. 

The voltage of the electrode is measured by defining 
a "measurement volume" near the back of the simula- 
tion cell — occupying the range — L/4 < a; < L/4, 
-L/4 < y < L/4, L/4: < z < 3L/8, where the origin 
{x, y, z) — (0, 0, 0) is located at the center of the electrode 
surface — inside of which the electric potential is mea- 
sured. After performing thermalization, the total elec- 
tric potential is measured at 500 equally-spaced points 
within the measurement volume after every 3(Ma -I- Mc) 
attempted moves. These measured values of potential are 
then averaged both temporally and spatially to produce 
a value for the voltage V of the electrode relative to the 
bulk. There was no noticeable systematic variation in 
electric potential across the measurement volume. The 
capacitance C{V — 0) is determined from the discrete 
derivative AQ/Ay at sufficiently small values of Q for 
which the relationship Q vs. V is linear. 

Our results are shown in Fig. [TJ For all Na^ stud- 
ied, the lowest value of T* in Fig. [1] is larger than the 
corresponding liquid-gas or liquid-solid coexistence tem- 
perature. For comparison, the triple point in the phase 
diagram of the RPM ionic liquid is located at No? — 0.5 
and T* = 0.025 while the gas-liquid critical point is at 
Na^ = 0.02 and T* = 0.05 :15^ We verified for each 
simulation that there was no phase separation within the 
simulation cell. 

In order to quantify the finite-size effects of our simu- 
lation cell, we examined the capacitance at zero voltage, 
C(0, T*), obtained from three simulation volumes of size 
L = 10a, 20a, and 30a. For Na^ = 0.4, C(0, T*) was seen 
to scale linearly with 1/L at all values of the tempera- 
ture that we examined (T* = 0.042, 0.072, and 0.14). In 
each case, the value of C{Q,T*) obtained by extrapola- 
tion to infinite system size was within 20% of the value 
of C(0,T*) corresponding to L = 20a. These results al- 
low us to conclude that the simulation cell with L = 20a 
provides a reasonable approximation of an infinite sys- 
tem. All MC results presented below correspond to this 
choice. 



IV. SEMI-QUANTITATIVE THEORY OF THE 
CAPACITANCE OF METALLIC ELECTRODES 

Our goal is to explain the large capacitance of the 
metal/ionic liquid interface at T* <C 1. As we empha- 
sized above, an ionic liquid next to a metallic electrode 
has a high degree of degeneracy because of the zero- 



energy process of ion pairs in the bulk dissociating and 
sticking to their images on the metal surface. As a re- 
sult, at low temperatures T* <C 1 effectively all ions in 
the system are either paired in the bulk or bound to 
their images on the metal surface [1^]. At zero applied 
voltage, equal numbers of positive and negative ions are 
bound to the metal surface. The area density ng of 
these ions can be estimated from the requirement that 
the chemical potential of pairs in the bulk be equal to 
the chemical potential of ions at the surfa ce, w hich gives 
ln{l/Na^) ~ 21n(l/noa2), so that uq ~ ^^N/a. 

As the voltage V of the electrode is increased from zero, 
some number of pairs in the system are separated so that 
the free counterion can come to neutralize the electrode 
surface. The corresponding density of these "excess ions" 
dn on the metal surface is related to the charge density 
a hy 6n — |cr|/e. If cr > 0, then Sn represents an excess 
of anions on the surface; if cr < the excess ions are 
cations. Naturally, excess ions condensed onto the metal 
surface will repel each other. Since each ion on the metal 
surface is separated by a distance a from its image charge 
in the metal, ions and their images constitute compact 
ion-image dipoles with dipole moment ea, and so excess 
ions repel each other via a dipole-dipole interaction 



u(6n) 



e^a^idn)^/^ 



2e 



(5) 



Excess ions at the metal surface are surrounded by no 
other ions per unit area, which effectively neutralize each 
other by forming no/2 bound pairs (see Fig. [5]). These 
no/2 bound pairs, along with bound pairs in the bulk, 
may serve to modify the effective dielectric constant for 
the interaction of excess ions. We comment on this pos- 
sibility later in this article. 

At low temperatures, the excess ions will seek to max- 
imize the distance from each other while maintaining a 
given density 5n, which results in the formation of a 
strongly-correlated liquid of excess ions reminiscent of 
a two-dimensional Wigner crystal (see Fig. [2]). The cor- 
responding total electrostatic energy per unit area of the 
system is 



U ~ a ■ 6n ■ u{5n) — aV, 



(6) 



where a is a numerical coefficient which describes the 
structure of the lattice of excess ions; for a triangular 
lattice, a ~ 4.4 jTZl]. The term —aV describes the work 
done by the voltage source. Fig. [2] shows a schematic 
depiction of the layer of ions bound to the metal surface. 
The voltage V which corresponds to a given charge 
density a can be found by the equilibrium condition 
dU/da = dU/d{eSn) ~ 0, which gives 
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The resulting capacitance per unit area C — da/dV — 
e[dV/d{6n)]''^ is 
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FIG. 2: A schematic depiction of the layer of ions bound to 
the metal surface at finite charge density a < 0. Cations 
(darkly-shaded spheres) and anions (lightly-shaded) are both 
bound to the metal surface by the strong attraction to their 
image charges. Excess cations (made even darker for em- 
phasis), which neutralize the electrode charge, minimize their 
repulsive energy by forming a Wigner crystal-like lattice. All 
other ions remain bound in neutral pairs. 



Substituting Eq. (O into Eq. 
terms of voltage: 



gives the capacitance in 
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This expression can be significantly larger than Ch when 
V is small. Physically, at such small voltages the excess 
ions are very sparse, and so their mutual repulsion goes 
to zero. In other words, at low voltages charging of the 
electrode is not limited by the accumulation of complete 
charged layers, but by the weak dipole-dipole repulsion 
between discrete, correlated ions and their image charges. 

Of course, the validity of Eqs. ([5]) - (O is limited to 
the range of voltage where there is a small fractional 
coverage of the metal surface by excess ions, na^ ^ 1. 
By Eq. ([7]), this corresponds to a dimensionless voltage 
V* = V/{e/ea) < 5a/4 « 5.5. At large enough voU- 
ages that na^ ~ 1, excess ions constitute a uniform layer 
of charge, and therefore the capacitance approaches Ch- 
At even larger voltages, the capacitance declines as com- 
plete layers of counterions accumulate next to the elec- 
trode and the double-layer becomes thicker. This leads 
to a mean-field capacitance C ex V~^^^ at large voltages, 
as described in Ref. f7i|. 

Formally, Eq. ^ diverges as the voltage goes to zero. 
Of course, this expression neglects entropic effects of the 
excess ions, which tend to destroy the lattice structure 
of excess dipoles on the metal surface. Such effects will 
truncate the low- voltage divergence of Eq. ^, result- 
ing in a finite capacitance at zero voltage. The value of 
this capacitance maximum can be estimated roughly by 
setting u(6nc) = ksT, solving for the corresponding con- 
centration Sric, and then plugging Suc into C{Sn) from 



Eq. 



This gives 
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(10) 
where A sa 0.6. In other words, the effective thickness 
d* = a{T*f/^/2A. At T* < 1, we find that d* < a. The 
corresponding voltage at which the capacitance plateaus 
{udd becomes equal to ksT) is V* = haT* /2 w IIT*. 

Fig. [3] shows the capacitance as a function of V* , as 
measured by our MC simulation, at density Na'^ = 0.4 
and at two values of the temperature T* . The inset shows 
the dimensionless charge density a* = aa^ /e of the elec- 
trode as a function of the voltage V* for the temperature 
T* = 0.042. The capacitance is determined by a numer- 
ical derivative of the a vs. V curve. Here, e/a^ is the 
maximal density for a square lattice of ions on the metal 
surface, so that a* can be interpreted as a filling factor 
of the first layer of ions. Note that the capacitance drops 
substantially even at low filling factor a* , so that the ca- 
pacitance is already reduced by a factor two at cr* = 0.5. 
This suggests that the decline in capacitance is not driven 
by the excluded volume effects emphasized in the theory 
of Ref. 0. 



action suggested in Eq. ([5]) is weaker at large ion density. 
One possible explanation is that at high densities ion 
pairs in the vicinity of two excess ions can polarize in the 
direction of the electric field, thereby producing an effec- 
tively larger dielectric constant. If we replace e in Eq. 
^ by an effective dielectric constant ee', then we find 
that A w 0.6(e')^/'^. The values of the constant A from 
above suggest that for bulk densities Na^ = 0.08,0.4, 
and 0.64, the value of e' is 1.0, 2.4, and 2.9, respectively. 
These values are consistent with our interpretation that 
the effective dielectric constant should increase with ion 
density, driving the capacitance upward. 

Based on our arguments from this section about the 
dependence of the capacitance on voltage and tempera- 
ture, we can hypothesize a general scaling relationship 
C{V*,T*) which reproduces Eqs. ® and ([TUll: 



C{V*,T* 



/3i 



p2T*)2 + (y*) 



211/6" 



(11) 



Here, /3i and (32 are numerical coefficients. Applying this 
fit to the curves shown in Fig. [3] gives a reasonably good 
fit with /32 « 8, suggesting that the capacitance plateaus 
at about V* = 8T*, as compared to the theoretically 
estimated value V:* = 11T*. 
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FIG. 3; The capacitance as a function of the dimensionless 
voltage V* — V/{e/ea) at two diflferent temperatures for a 
system with ion density Na^ = 0.4. The inset shows a plot of 
the dimensionless charge density a* — aa^/e as a function 
of the voltage V* measured by the MC simulation at the 
temperature T* = 0.042. For e = 3 and o = 1 nm, V* ^ 2 
corresponds to 0.96 Volts. 



The prediction of Eq. ([T0| provides a good fit to the ca- 
pacitance measured by our MC simulation at low ion den- 
sity, as shown in Fig. [1] However, while the dependence 
C ex (T*)~^/^ remains accurate for all densities, the con- 
stant A apparently depends on the ion density, taking 
the values A = 0.6, 1.1, and 1.2 for Na^ = 0.08,0.4, and 
0.64, respectively. This increase is also reflected in Fig. 
[3l where the capacitance at finite voltage is somewhat 
larger than predicted by Eq. ©, consistent with the in- 
crease in the constant A. This larger capacitance at high 
densities is perhaps an indication that the dipole inter- 



V. ELECTRODE MATERIAL: FROM PERFECT 
TO POOR METAL 

So far we have assumed that the electrode is a perfect 
metal, or in other words, that the screening radius b of 
the metal is much smaller than the ion diameter a. This 
assumption is justified for ionic liquids with large ions 
and electrodes made from a good metal. Experiments 
on such systems have indeed reported large capacitance 
that declines with absolute value of voltage [l^ (the C- 
V curve is "bell-shaped"). However, for smaller ions and 
for electrodes made from semi-metals, such as graphite 
or glassy carbon, experimental values of C(0) are smaller 
and the C{V) curves are "camel-shaped" , i. e. the capac- 
itance grows parabolically near V = |18l - [20J . 

In order to interpret this difference qualitatively, let 
us recall that in such poor metals the density of states 
at the Fermi level is relatively small and the screening 
radius b of the metal may become comparable to a/2. As 
a result, the image potential may change. When b < a/2 
one can think that the electric field produced by ions at 
the metal surface is relatively weak and slowly-varying. 
In such a case the screening charge of the metal is ef- 
fectively situated at the distance b away from the metal 
surface, i.e. at z ~ —b. Therefore, the reflection plane 
for the image charge is at z = —6, so that an ion at dis- 
tance z from the surface experiences a smaller attraction 
~e'^/4e{z + b) to the surface, rather than the standard 
—e'^/Aez for a perfect metal. At the distance of closest 
approach z = a/2, the ion-to-surface attraction energy 
becomes —e'^/2e{a + 2b). This leads to a finite energy 
cost Eq for dissociating a bulk ion pair and bringing it to 



the metal surface, given by i^o = e^ / e[a~^ — {a+2b)~^]~^ . 
Thus, a finite voltage is necessary to break pairs in the 
bulk and obtain free ions which can provide screening. 
This means that, for an electrode with finite screening 
radius b, the bell-shaped C{V) curve splits into two peaks 
located at V^ = zLEq/c, thereby becoming camel-shaped, 
in agreement with the above-mentioned data. 

The capacitance C{V) at F > Eo/e can be estimated 
with the help of the theory in section IV leading up to Eq. 
©. In this case, however, the voltage V in Eq. © should 
be replaced by y — Eq/c and the dipole arm a should be 
replaced by the longer arm a + 26. These substitutions 
give 
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e{a + 2b)iV-Eo/e) 



1/3 



26 



Ch (12) 



for V > Eo/e. Since the dipole-dipole repulsion is 
much stronger due to the longer dipole arm, C{V) is 
substantially smaller and reaches its geometrical limit 
C_f/(6) — Ch ■ a/ {a + 26) at a smaller voltage V — Eq/c — 
5ae/4e{a + 26), or V* - Eo/ie^/ea) w 5.5/(1 -I- 26/a). 
Starting from this voltage the capacitance saturates at 
the level of Cff (6). 

Only at even larger voltage T^* — £'o/(e^/ea) > [5.5a^-|- 
8TTb{a + b)]/[a{a+2b)] do counterions comprise a full layer 
at the surface, after which the theory of multi-layer ar- 
rangement of ions [71 becomes applicable. This same be- 
havior for C{V) is expected in the case where a good 
metal is covered by a thin insulating layer, for example, 
its own oxide. 

In order to verify this theory we repeated our MC cal- 
culations for T* — 0.04 and Na^ = 0.5 using a relatively 
large 6 = a/2, which is at the limit of applicability for 
linear screening by the electrode surface. For simplicity, 
we have also assumed that the metal ion lattice has the 
same dielectric constant as our ionic liquid. Results are 
shown on Fig. lU plotted as a function of V* and a* . As 
expected, the peak at F = is split into a camel-like 
structure (we show only the positive half of the symmet- 
ric C-V curve). The characteristic dimensionless voltage 
of the peak is V* ~ 0.5, in agreement with the above esti- 
mate for Eq . Note that the capacitance maximum occurs 
at cr* = 0.1 and is apparently not related to excluded vol- 
ume effects among counterions. As predicted by Eq. (|12|) . 
the peak capacitance C{V*) is approximately 2.5 times 
smaller than at 6 = (recall that CH{b) in Fig.|3]is twice 
smaller than Ch in Fig. [3]). 

Figs. 1-3 clearly show that, for both a perfect metal 
(6 = 0) and a semi-metal with 6 — a/2, the capacitance 
C{V) can be 2-3 times larger than the geometrical mean- 
field capacitance Cnib) at moderately low temperature 
and voltage. One can interpret this fact by saying that 
the geometrical capacitance is in series with a negative 
capacitance from the ionic liquid, Cu, so that C^^ = 
Cff(6)-i +C~i^ < Cniby^- Multiplying this equation 
by e/An we arrive at d* = a/2 + b + r^, where r^ is 
the screening radius of the ionic liquid. As we saw, at 
low temperatures and voltages this definition leads to 




FIG. 4: The ratio of the capacitance of the metal/ionic liquid 
interface C to the geometrical capacitance Cnib), plotted as a 
function of dimensionless voltage V* (bottom axis) and charge 
density a* (top axis) for a system with ion density Na^ = 
0.4 and metal screening radius b = a/2. The capacitance 
is determined by a numerical derivative of the a vs. V data 
obtained from a MC simulation at temperature T* — 0.042. 
For £ = 3 and o = 1 nm, V* — 4 corresponds to 1.92 Volts. 



negative Tj. At high temperatures T* > 1, the screening 
radius is positive and given by Eq. (U}: r^ = toh- Then, 
as explained in the introduction, d* ~ r^. 



VI. ANALOGY WITH SEMICONDUCTOR 
DEVICES 

Capacitance smaller than the geometrical one, or in 
other words negative screening radius r^, is well-known 
in semiconductor physics for capacitors made of a metal, 
an insulator of width d, and a semiconductor contain- 
ing a clean two-dimensional electron gas (2DEG) with 
two-dimensional density n. This can be, for example, a 
Si MOSFET or a gated GaAs-GaAlAs heterostructure. 
In the limit of low density n, a 2DEG is a classical sys- 
tem whose physics is dominated by the Coulomb interac- 
tion between electrons, leading to a Wigner crystal-like 
strongly-correlated liquid state. This state was shown 
to have a negative thermodynamic density of states and 
a negative screening radius rg — — 0.116n~^/^. When 
77,-1/2 ^ ^^ |.]^jg gcpegjiing radius produces a small nega- 
tive correction to the geometric capacitance width d, so 
that HH-il 



d-0.116n 



-1/2 



(n-1/2 « ay 



(13) 



The question of what happens to d* when n^^/^ ^ d 
has never been addressed. Now we understand that in 
this limit one should think about the 2DEG and its im- 
age charges as a gas of electron-image dipoles oriented 
along the z axis with dipole arms 2d (we consider the 
electrons, unlike ions, to be point-like particles so that 



for this case a — 0). Repeating the calculations leading 
to Eq. © , again using the point-dipole approximation for 
the dipole-dipole repulsion, leads at low T to an anoma- 
lously large capacitance C{V) corresponding to the effec- 
tive capacitor thickness 



d* = 2.68d\ 



(n-i/2 » d). 



(14) 



In order to study the crossover between Eq. (fT3| and 
Eq. (fT4|) . we calculated a universal low temperature func- 
tion d* {d, n) through an exact numerical calculation of 
the interaction energy of an infinite lattice of dipoles with 
dipole arm 2d and areal density n, abandoning the point 
dipole approximation. The result can be written as 



d-f{dn^'^) 



(15) 



where f{x) is a dimensionless function shown in Fig. [5] 
In the case of a low-temperature ionic liquid with diam- 
eter small enough that na^ <C 1, one can use this result 
by replacing d with d + a/2. Then Eq. p^ describes 
the shape of the crossover of C{V) to the geometrical 
capacitance C/f (6). 




FIG. 5: The scaling function f{x) which determines the ef- 
fective thickness d* of a capacitor composed of a 2DEG with 
area density n separated from a perfect metal surface by an 
insulating layer of thickness d, as defined by Eq. (jlSp . The 
left side of the plot corresponds to a very sparse 2DEG, where 
the electrons can be thought to form an array of discrete ion- 
image dipoles and d* is described by Eq. (IT4)l . The right side 
corresponds to a relatively dense packing of electrons, where 
the electrons approach a uniform layer of charge and d* is 
described by Eq. ([13} . 



We note that the predictions of this section might be 
verifiable in an extremely clean GaAs-GaAlAs hetero- 
j unction with an area density of holes lower than n = 10^ 
cm~^ [2^, assuming that one is able to make gates even 
closer than the current limit of 250 nm without an in- 
crease in disorder. Our theory might also be verified in 
a capacitor composed of a 2DEG on top of a thin film of 
liquid Helium covering a metal gate 



VII. CONCLUSION 



To summarize, this paper is concerned with the re- 
stricted primitive model of an ionic liquid at a metal 
interface. Within this model, we obtain capacitance at 
zero voltage as large as 3C_f/- We also find that for a good 
metallic electrode at small voltage, C decreases with T as 
r^i/3_ When the temperature is fixed and is relatively 
small, C{V) decreases as 1/V^^^ (the C-V curve is "bell- 
shaped"). On the other hand, when the electrode is made 
from a semi-metal the C-V curve is "camel-shaped", 
meaning that the capacitance first grows with V and then 
goes through a maximum and decays as 1/V^''^ . We in- 
terpret these results with the help of a semi-quantitative 
analytical theory based on the weak repulsion between 
ion-image dipoles, and we confirm our results with a MC 
simulation. Our conclusions are in qualitative agreement 
with experimental data. 
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